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Abstract 



An alternate formalism is developed to determine the energy eigenvalues of quantum 
mechanical systems, confined within a rigid impenetrable spherical box of radius r , in the 
framework of Wentzel-Kramers-Brillouin (WKB) approximation. Instead of considering 
the Langer correction for the centrifugal term, the approach adopted here is that of Hainz 
& Grabert : The centrifugal term is expanded perturbatively (in powers of h), decompos- 
ing it into 2 terms — the classical centrifugal potential and a quantum correction. Hainz 
and Grabert found that this method reproduced the exact energies of the hydrogen atom, 
to the first order in h, with all higher order corrections vanishing. In the present study, 
this formalism is extended to the case of radial potentials under hard wall confinement, to 
check whether the same argument holds good for such confined systems as well. As expicit 
examples, 3 widely known potentials are studied, which are of considerable importance in 
the theoretical treatment of various atomic phenomena involving atomic transitions, viz., 
the 3-dimensional Harmonic oscillator, the hydrogen atom and the Hulthen potential. 



key words : WKB approximation, 3-dimensional spatial confinement, radial potentials, 
perturbative expansion, centrifugal term, harmonic oscillator, hydrogen atom, Hulthen 
potential 
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1. Introduction 



Spatial confinement of electrons in artificial nanostructures, on a scale comparable to 
their de Broglie wavelength, is a much talked about subject for the past decade or so 
[1-7]. Such artificial atoms as they are called because of their quantized energies, undergo 
radical changes in terms of both physical and chemical properties, because of their ex- 
tremely small spatial dimensions, making them very useful in the study of atomic and 
molecular phenomena [8]. However, the problem with these so-called quantum wells, 
quantum wires and quantum dots is that their exact analytical treatment is not possible 
in most of the cases. Consequently, various approximation methods — the variational ap- 
proach, the shifted 1/N expansion technique, the modified Airy function (MAF) method, 
the supersymmetric version of the same (SMAF), the Wentzel-Kramers-Brillouin (WKB) 
approximation, its supersymmetric version (SWKB), etc. — come into the picture [9-18]. 
The other option is to go for a numerical solution. Of the several approximation meth- 
ods, the semiclassical WKB approximation technique is a very effective tool, yielding 
fairly accurate results in various quantum mechanical problems, where exact solutions are 
unknown or difficult to find out. This method gives a good estimate of both the energy 
eigenvalues as well as eigenfunctions, with the exception of the region near the classical 
turning points. Though this approach works well for one-dimensional problems, practical 
use shows that the standard leading order WKB approximation always reproduces the 
exact spectrum for the solvable spherically symmetric potentials V(r) if the centrifugal 
term 

Voir-) = 1(1 + l)h 2 /2r 2 
is replaced by the Langer correction term [19] 

V L (r) = (/ + l/2)V/2r 2 

This modification also justifies the WKB expansion of singular potentials like that of the 
Coulomb potential, near the origin. However, some authors have attempted to get rid 
of this Langer modification (LM), based on non-linear transformations [20] and super- 
symmetry [21]. Though their approach yielded results which are superior to those with 
LM, Coulomb type problems did not fare well. Hainz and Grabert [22] challenged this 
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common belief and put forward a new method to deal with centrifugal terms in the WKB 
approximation. 

Since the semiclassical WKB approximation proceeds as a perturbation in powers of 
h, it was argued in ref. [22] that within this expansion, the centrifugal term can be 
decomposed as 

l(l + l)h* _ Ll | hL 
2mr 2 2mr 2 2mr 2 
with Lq = hi. The first term is the classical centrifugal term, while the second term is a 
quantum correction. Thus the quantum correction can be treated as a perturbation and 
expanded accordingly. Proceeding along these lines, Hainz and Grabert [22] found that 
the semiclassical energy eigenvalues for the hydrogen atom turned out to be exact to the 
first order in h, with all higher order corrections vanishing. The aim of the present work 
is to check whether the same argument holds good for confined systems as well. The 
formalism of ref. [22] is extended to the case of radial potentials confined within rigid 
impenetrable spheres of radius r . As expicit examples, 3 widely known cases are studied, 
viz., 

(i) the 3-dimensional Harmonic oscillator (HO) 

(ii) the hydrogen atom 

(iii) the Hulthen potential 

These potentials are of considerable importance in theoretical treatment of various atomic 
phenomena involving atomic transitions. It is observed that this formalism gives a better 
estimate of the energy spectrum than the case with Langer modification, even in case of 
spatially confined systems. 

The organisation of the paper is as follows. In Section II, the WKB formalism is ex- 
tended to include quantum mechanical systems confined radially, expanding the centrifu- 
gal term perturbatively in powers of h. In Section III, the approach is applied explicitly 
to 3 physically relevant cases, viz., the 3-dimensional harmonic oscillator, the hydrogen 
atom, and the Hulthen potential. Section IV is kept for conclusions and discussions. 
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2. Theory 



The starting point of the study is the three-dimensional Schrodinger equation for a radial 
potential V(r) 

l{l + l)h 21 



d 2 ,. , 2m 

dV 2 m + ?T 2 



E 



V(r) 



2mr 2 







(2) 



where 1(1 + l)h 2 represents the eigen values of the square of the angular momentum 
operator L 2 and m is the mass of the particle. It is worth noting here that the WKB 
approximation can be applied only when the de Broglie wavelength A = h/p (h — 2irh) 
is changing slowly. With the help of (1) the radial Schrodinger equation (2) can be cast 
in the form 



dr 2 



2m 



E 



where 



V ef f(r) 
V eff (r) = V(r) + 



h 



2mr 2 



ip(r) = 



T 2 

2mr 2 



(3) 



(4) 



In order for the physical system to have a stable bound state (discrete spectrum) it must 
have two classical turning points r\ and r 2 . This gives rise to 3 regions given by 



region I 
region II 
region III 
where 



< r < n 
r\ < r < r 2 
r > r 2 



V 1 > E 

E > V\ 

V 1 > E, 

V 1 (r) = V eff (r) + h 



2mr z 



If one defines T(r) and re(r), with (k 2 = —T 2 ) by 



then (3) reduces to 



dr 2 



- K 2 ( r ) = 



in regions I and III 



(5) 



(6) 



(7) 



(8) 
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{* + iV)L=o 



in region II 



(9) 



Expanding in powers of H, and keeping terms to the first order in h, one can write T(r) 
and k{t) as 

r(r) ~ r (r) - 



/c(r) ~ /«o( r ) + 



2hT r 2 



where 



2m 

Ur) = J W {E-V eff } 



so that 



^{r) = f-^{V eff -E} 



1 + 



i i 



4/U> 2 



The conventional WKB ansatz is assumed for the wave function 



■0(r) = exp 



^E(-^) fe 5 fc (r) 



Substituting 



Vk(r) = 



dS k (r) 
dr 



and expanding them in powers of H, one obtains the set of relations 

yo = ±y/2m(E-V eff (r)) 



1 ( , , .L 



2y 



2m-2 



V2r 



~2y~ \f m + y ' 2m ~ l + 2 ^ V2m ~ k Vk 



k=i 



2m- 1 



Z/2m+l = ~— W 2 m + 2 V^m+l-k Vk 



(10) 
(11) 

(12) 
(13) 

(14) 
(15) 



(16) 
(17) 

(18) 
(19) 

(20) 
(21) 



k=i 
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In the above relationships, prime denotes differentiation with respect to r. Thus y turns 
out to be the classical momentum, and 
y = ±iK h in regions I and III, 
yo = T h in region II. 

So, the wave function is a linear combination of the form 



where 



■^Jdr y ± (r) 



(22) 



y(r) = E(-ih) k y k (r) (23) 

with derivatives (to the first order in %) 

dip ( i 



dr + (24) 



d 2 tp S ' 2 i 



+ - (2S 'S 1 ' + S ' ') + ^" W> (25) 



dr 2 \ h 2 ' h 

This gives the complete solution to the Schrodinger equation ( to the first order in ft, ) as 

ip = — U exp\- [yodr- 1 -^- f 1 (26) 
v^o [h J 2 J y r 2 J 

Since the radial wave function must vanish at r = 0, the only allowed solution in region 
I is 

fa = ^=exp I - / n dr - -| / \ (27) 

v /fi^ [ Jr 2Tl Jr K r 2 J 

Now we shall make use of the conventional connection formulae for WKB approximation 
at the turning point r±, [23] viz., 

exp ^— J n(r) dr^j = sin^J T(r) dr + — j (28) 



J4r) V Jr ) fifi 

l / n , , . \ l 

exp 



(f K{r)dr ) = jm cos {i T(r)dr + j) (29) 



'/c(r) Wr / y/r(r) 

Expanding T(r) and k(t) in terms of To(r) and Ko(r), and keeping terms upto the first 
order in ft only, the connection formulae can be cast in the form 



^m exp (iy (T)dT + Tn[ 



(30) 



ft (r)r 2 / 

*HjC™*-5jCw + 5} (31) 

The solutions to the Schrodinger equation in regions II and III are obtained by matching 
the WKB solutions on either side of the turning points r± and r 2 , with the help of the 
connection formulae (30) and (31). Thus the solution in region II comes out to be 

To obtain the solution in region III, ipu is written as 

= £ To(r)rfr - I? T (r)dr - % ^ + J? ^ + f 
= (f + £ To(r)rfr - § £ ^) - (/- r (r)dr - § ^ + 
= (f + 6l)-S 



(33) 



where 



and 



rr2 „ , x , £n f 2 dr 7T 



B = i r «M*-^/ f^ + I < 35 > 

Substituting (33) into (32) and simplifying, 

2A 

ip n = (cos 9 cosB + sin 9 sin B) (36) 

\/ r oW 

Using the connection formula at the turning point r 2 , the solution in region III is obtained 

as 

2A A 
ipjjj = , cos 9 exp(cr(r)) + . sin 9 exp(—cr(r)) (37) 

\/ K o( r ) \JK {r) 
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where 

*> = jO (r) *-5s££ (38) 

Now, the WKB quantization rule is obtained by the constraint ip(r = r ) = 0, where 
r is the radius of the confining spherical box. 2 cases arise depending on the size of 
confinement, i.e. 

i) there is a single turning point within the box ( T\ < r < r 2 ) 

ii) both the turning points are within the box (r > r 2 ) 
thus yielding 2 different quantization rules. 



1) Extremely small confinement : 

The size of the box is so small that it admits only a single turning point, i.e., r\ < r < r 2 
This modifies the boundary condition to : ipn(ro) = 0. 
Hence the WKB quantization condition reads 

Ai - A 2 = (n+^), ra = 0,l,2, (39) 



with 

X, = / T (r)dr (40) 
L r r ° dr 



A2 = fnL f^F (41) 



2) The confining box is not so small : 

The size of the box is such that both the classical turning points lie within it, i.e., r > r 2 

Hence, in this case the solution of the Schrodinger equation must obey the boundary 
condition : ipni(r ) = 0, 
yielding the WKB quantization rule 

2 cos 9. exp(a(ro)) + sin 9. exp (— o"(r )) = (42) 

These quantization rules (39) and (42) enable one to determine the energy spectrum of 
any spatially confined, radial potential, in the framework of WKB approximation. 
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3. Calculations 



In this section, the WKB rules so developed are used to determine the energy spectrum 
of 3 explicit potentials 

(i) the 3-dimensional Harmonic oscillator 

(ii) the hydrogen atom 

(iii) the Hulthen potential 

each confined within an impenetrable spherical box of radius r . All the 3 potentials are of 
tremendous importance in a variety of physical problems, and have been studied widely 
Units used are h — m — 1 so that L = I. 

3-dimensional Harmonic Oscillator 

V(r) = j (43) 
Hence the relationship E — V e ff = gives the classical turning points at 

n = { E -Ve 2 -1 2 } 1/2 (44) 

r 2 = [E + VW^P} 112 (45) 

Proceeding along the formalism developed above, the energy levels of the enclosed 3- 
dim. harmonic oscillator are computed for various values of the confining parameter r , 
with the help of the mathematical relationships in ref. [24]. The results are presented 
in Table I, comparing the energies so obtained (with no modification of the centrifugal 
term) denoted by E, with those from the conventional WKB quantization rules for 3- 
dimensional confinement (with Langer modification), E(WKB) [18], the direct variational 
method, E(var) [9], and exact numerical values, E(exact) [9] for the enclosed 3-dim. 
harmonic oscillator. 

Table I : Enclosed 3-dim. Harmonic Oscillator 

(n — n r + I + 1 : n r — 0, 1 — 1) 
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E 


E(WKB) 


E(var) 


E(exact) 


1.0 


10.2876 


10.2643 


10.3188 


10.2822 


1.5 


4.9068 


4.9084 


4.9169 


4.9036 


2.0 


3.3081* 


3.2490 


3.2514 


3.2469 


2.5 


2.6835 


2.7079* 


2.6901 


2.6881 


3.0 


2.5313 


2.5310 


2.5337 


2.5313 


4.0 


2.5001 


2.5001 


2.5015 


2.5001 


5.0 


2.5000 


2.5000 


2.5012 


2.5000 



* In these cases the size of the rigid spherical box is such that the wall is close to the 
turning point, where the WKB approximation is not expected to give good results. 



Hydrogen Atom 

The well-known Coulomb potential 

V(r) = - - (46) 
r 

is known to possess negative energies. However, spatial confinement alters this scenario. 
For extremely small confinement, the system is no longer bound. E > 0, and E—V e ff = 
gives a single turning point at 

n = (47) 

For all practical purposes, r t is very small, and the eigenenergies of the enclosed system 
are obtained from the relationship (39) 

r° i I r° dr / 3\ 

/ J2 (E - Veff) dr - - / . = (n + ^)tt (48) 

Jrt V 1 //J 2 Jr t p (E-V eff ) r 2 V 4^ K J 

with 

Veff = -\ + ^ (49) 

However, for bound energies (E < 0), 2 cases may arise, depending on the size of the 
confining box : 

Either there is only one turning point inside the box (ri < ro < r2), or confining wall 
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encloses both the turning points (r > r 2 ). Let E = — e and V e ff = —v e ff. The 
expression E — V e ff = gives the roots at 



The energy spectrum of the boxed-in hydrogen atom is determined with the help of the 
formalism developed above. The results are computed and presented in tabular form, 
for the 2p (Table II) and 3d (Table III) states. The energy eigenvalues calculated in 
this study, E, are compared with those obtained by the conventional (with LM) confined 
WKB approximation developed earlier E(WKB) [18], the direct variational method of 
Marin and Cruz, E(var) [9], Varshni's modification of Marin-Cruz approach, E(Varshni) 
[25], and the exact numerical values, E(exact) [25,26]. 




(50) 




(51) 
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Table II : Enclosed Hydrogen Atom - 2p state : (n — n r + 1 + 1 : n r — 0, 1 — 1) 



r 


E 


E(WKB) 


E(var) 


E(Varshni) 


E (exact) 


0.6 


49.8448 


49.3997 


50.401 


49.935 


49.874 


0.8 


26.9179 


26.5586 


27.155 


26.910 


26.879 


1.0 


16.5063 


16.2590 


16.611 


16.464 


16.446 


1.2 


10.8828 


10.7653 


10.999 


10.905 


10.893 


1.4 


7.6209 


7.4379 


7.6857 


7.6214 


7.6138 


1.6 


5.5112 


5.3928 


5.5801 


5.5347 


5.5295 


1.8 


4.1512* 


4.0693 


4.1675 


4.1345 


4.1308 


2.0 


3.1513 


3.1010 


3.1791 


3.1547 


3.1520 


2.2 


2.4469 


2.4013 


2.4641 


2.4458 


2.4438 


2.4 


1.9224 


1.8815 


1.9326 


1.9187 


1.9173 


2.8 


1.2129 


1.1807 


1.2157 


1.2075 


1.2068 


3.0 


0.9684 


0.9420 


0.9694 


0.9631 


0.9625 


3.5 


0.5466 


0.5371 


0.5459 


0.5427 


0.5424 


4.0 


0.2894 


0.2771 


0.2888 


0.2872 


0.2871 


5.0 


0.0154 


0.0135 


0.0155 


0.0152 


0.0152 


7.0 


-0.1687* 


-0.1666 


-0.1748 


-0.1748 


-0.1749 


10.0 


-0.2269 


-0.2256 


-0.2369 




-0.2377 


14.0 


-0.2487 


-0.2484 


-0.2484 




-0.2491 



Table III : Enclosed Hydrogen Atom - 3d state : (n = n r + / + 1 : n r = 0, 1 = 2) 





E 


E(WKB) 


E(var) 


E(Varshni) 


E(exact) 


1.0 


29.8203 


29.7306 


30.234 


29.979 


29.935 


1.5 


12.5321 


12.4895 


12.692 


12.587 


12.570 


2.0 


6.6415 


6.6064 


6.7182 


6.6640 


6.6550 


2.5 


3.9882 


3.9658 


4.0288 


3.9970 


3.9920 


3.0 


2.5863 


2.5593 


2.6088 


2.5887 


2.5856 


4.0 


1.2467 


1.2379 


1.2532 


1.2440 


1.2427 


5.0 


0.6581 


0.6489 


0.6634 


0.6588 


0.6582 


6.0 


0.3617 


0.3550 


0.3634 


0.3609 


0.3607 


7.0 


0.1926 


0.1890 


0.1945 


0.1933 


0.1932 


8.0 


0.0919 


0.0897 


0.0928 


0.0922 


0.0921 


10.0 


-0.0140 


-0.0156 


-0.0141 


-0.0142 


-0.0142 


12.0 


-0.0625 


-0.0626 


-0.0625 


-0.0625 


-0.0625 


14.0 


-0.0862 


-0.0860 


-0.0862 




-0.0862 


16.0 


-0.0939* 


-0.0928 


-0.0982 




-0.0984 


20.0 


-0.1079 


-0.1077 


-0.1076 




-1.1079 
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Confined Hulthen potential : 

Screened Coulomb potentials are of tremendous importance in atomic phenomena. The 
particular case studied here is the confined Hulthen potential, given by 



V (r) = V -^-^ (52) 
where Vq = —Z5, with Z the atomic number and S the screening parameter. Taking 

z = i, 

V ef f(r) = V(r) + ^ (53) 

Once again the energy eigenvalues are computed for different values of the confining 
radius r and screening parameter 5, and the results presented in tabular form, for ease 
of comparison with those obtained from other approximation methods, viz., the shifted 
1/N expansion method E(l/N) [15], and exact numerical energies E(exact) [15]. 

Table IV : Confined Hulthen potential : (5 = 0.1, n = n r + I + 1 ) 





state 


n r 


I 


E 


E (exact) 


E(l/N) 


6 


2p 





1 


-0.00782 


-0.00865 


-0.00294 


7 


2p 





1 


-0.03976 


-0.04069 


-0.03324 


8 


2p 





1 


-0.05510 


-0.05783 


-0.05293 


9 


2p 





1 


-0.06612 


-0.06728 


-0.06389 


10 


2p 





1 


-0.07196 


-0.07257 


-0.07008 


25 


2p 





1 


-0.07921 


-0.07918 


-0.07920 




3p 


1 


1 


-0.01384 


-0.01475 


-0.01295 




3d 





2 


-0.01381 


-0.01390 


-0.01332 


50 


2p 





1 


-0.07920 


-0.07918 


-0.07920 




3p 


1 


1 


-0.01598 


-0.01605 


-0.01578 




3d 





2 


-0.01450 


-0.01448 


-0.01450 


Table V : 


Confined Hulthen potential : (5 = I 


r 


state 


n r 


/ 


E 


E (exact) 


E(l/N) 


8 


2p 





1 


-0.01607 


-0.01731 


-0.01242 


9 


2p 





1 


-0.02612 


-0.02749 


-0.02428 


10 


2p 





1 


-0.03389 


-0.03339 


-0.03118 


25 


2p 





1 


-0.04192 


-0.04188 


-0.04199 


50 


2p 





1 


-0.04191 


-0.04189 


-0.04196 



.2,n = n r + I + 1 ) 
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4. Discussions and Conclusions 



In the present study the WKB approximation technique is used to derive an alternate 
formalism for quantum systems with radial potentials, confined within a rigid spherical 
box of radius r . Instead of the conventional Langer modification, in this approach the 
centrifugal term is decomposed perturbatively (in powers of H) into 2 terms — , the clas- 
sical centrifugal potential and a quantum correction, following the analysis of Hainz and 
Grabert [22]. The unique advantage of this approach is that it requires no modification 
of the centrifugal term in the WKB expansion when applied to radial potentials. More- 
over, the quantization rules follow naturally from the WKB connection formulae, and the 
calculations are straightforward, though somewhat lengthy. 

As a testing ground, the analysis is applied explicitly to 3 widely studied confined 
systems, each of considerable importance in atomic phenomena, viz., the 3-dimensional 
Harmonic oscillator, the hydrogen atom and the Hulthen potential. Each system is con- 
fined within a rigid spherical box of radius r . The spatial confinement imposes the 
additional boundary condition ip(r) = at r = r on the radial wave function. 
This criterion, alongwith the WKB connection formula, gives the quantization rules for 
estimating the energy eigenvalues, E. For each case the results are computed and pre- 
sented in tabular form, in Tables I - V, for ease of comparison with those obtained from 
other approximation methods, viz., the conventional WKB method for confined systems 
(using the Langer modification of the centrifugal term [18], variational results of Marin 
and Cruz [9], modified form of the same as given by Varshni [25] (in case of confined 
Hydrogen atom ), shifted 1/N expansion method [15] (in case of the confined Hulthen 
potential), and exact numerical results. 

It is easy to observe from the Tables I to V, that the present formalism works quite 
well for all the 3 cases. The energy values are better than the conventional WKB ener- 
gies (with Langer modification), as well as the shifted 1/N expansion results, justifying 
the perturbative expansion of the centrifugal term even for potentials under hard-wall 
confinement. In some cases these energies are even better than the variational results of 
Marin and Cruz [9] . This is true for most of the confining radii, except when the size of 
the box is close to the turning point (marked by * in the tables). This is expected as the 
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WKB approximation is not valid close to the turning points. It may be worth mentioning 
here that Hainz and Grabert version of the WKB method addresses primarily problems 
with singular potentials. Though the harmonic oscillator does not fall in this class strictly, 
neverthelass, the centrifugal term introduces a singularity at the origin for non-zero I. So 
the author found it interesting to check the behaviour of the confined harmonic oscillator 
under such an expansion. 

Another point worth examining is the effect of the higher order terms. It was shown in 
ref. [22] that in contrast to Langer modification, the higher than first order terms gave 
a vanishing contribution to the estimate for energy. However, this fact does not hold for 
confined potentials as is evident below : 
Expanding in powers of H, one obtains 

fro rro rr [ jn rr £2 . rQ r3 

/ T(r)dr= T (r)dr- —±-dr-l 9 ° - dr - / (54) 

J ri Jri Jri 41 IL o' Jri 

Proceeding along analogous lines as above it was found that the eigenvalues obtained in 
the lowest order became worse when higher order corrections were evaluated. 

To conclude, this present analysis of perturbative decomposition of the centrifugal term 
into 2 parts — a classical potential and a quantum correction — plays a vital role in 
improving the WKB quantisation rule in the first order, thus yielding better estimates of 
the energy eigenvalues. Hence this formalism may be useful in determining the energy 
spectrum of any 3-dimensional radially confined problem. 
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